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^ Abstract 

^ Optical chaos is a topic of current research characterized by high-dimensional nonlinearity which is attributed to the delay- 

induced dynamics, high bandwidth and easy modular implementation of optical feedback. In light of these facts, which 
adds enough confusion and diffusion properties for secure communications, we explore the synchronization phenomena in 
spatiotemporal semiconductor laser systems. The novel system is used in a two-phase colored image encryption process. The 
high-dimensional chaotic attractor generated by the system produces a completely randomized chaotic time series, which is ideal 
in the secure encoding of messages. The scheme thus illustrated is a two-phase encryption method, which provides sufficiently 
high confusion and diffusion properties of chaotic cryptosystem employed with unique data sets of processed chaotic sequences. 
' ^ ' In this novel method of cryptography, the chaotic phase masks are represented as images using the chaotic sequences as the 
elements of the image. The scheme drastically permutes the positions of the picture elements. The next additional layer of 
' |"i security further alters the statistical information of the original image to a great extent along the three-color planes. The 

O I intermediate results during encryption demonstrate the infeasibility for an unauthorized user to decipher the cipher image. 

O Exhaustive statistical tests conducted validate that the scheme is robust against noise and resistant to common attacks due to 

^ the double shield of encryption and the infinite dimensionality of the relevant system of partial differential equations. 
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I. INTRODUCTION 



The immense proliferation of internet and the advancement in telecommunication has led to vulnerability problems 
in securely communicating the information. Confidentiality of data and information are thus an area of concern 
for many researchers. An efficient and faster new breed of encryption technique operating in the physical layer 
of the transmission system is lately creating a lot of buzz. Quantum Cryptography and Chaos Cryptography are 
subtle deviations from the traditional methods of encryption techniques which are lately among the top research 
fiurry. The first primarily focuses on quantum key distribution. The second method deals with the scrambling and 
transfer of encoded messages at a very high encryption data rate up to several tens of Gb/s, unlike other software 
based traditional schemes. This method of chaos cryptography is the outcome of the phenomenal work of chaos 
synchronization demonstrated by Pecora and Carroll [1]. This bifurcation from the traditional methods of encryption 
realized with the synchronization of two coupled chaotic trajectories opened a whole new dimension for cryptographic 
applications [2-7]. Its principle of operation relies primarily on the fact that the chaotic nonlinear oscillator plays 
the role of a broadband signal generator and the chaotic waveform can mask the information. The most distinct 
feature of chaotic dynamics is that they possess ergodicity [8], sensitivity to initial condition, randomness and mixing. 
These are the key properties which are exploited in communication theory for secure transmission of data. The 
experiment conducted by Cuomo and Oppenheim [9] demonstrated the feasibility of such a communication scheme 
suggesting that the unpredictability of chaotic oscillations of synchronized systems can be used effectively to encrypt 
information-bearing signals, while their deterministic nature can be used in decryption. In communication theory, a 
waveform which is coupled with another can be formed as a carrier over the communication channel. Here, due to 
chaos synchronization, the receiver end can retrieve the embedded information. 

On the other hand, in order to design an efficient cryptographic system, it is imperative to select a system which is 
of high dimension and remains chaotic in a continuous range of the parameter space [10]. Moreover, the encryption 
methodology used to develop an efficient cryptic system should incorporate properties like a huge key space [11], 
which should be very sensitive to any kind of perturbations and immune from statistical attacks. Chaos based 
cryptography renders the cipher deterministically disordered having a strong dependence on even minimal variations 
of initial conditions and parameter values, and makes them resistant to most of the attacks. 

In apropos to the above discussion, another secure method of communication can be realized at the physical layer 
itself by using chaotic carriers obtained from semiconductor lasers to encode an information. Reference [12] elicits 
that these nonlinear systems on synchronization have tremendous potential to achieve enhanced security where they 
play the role of emitter and receiver system when subjected to a high injection strength [13]. It is well known that, 
higher the complexity of the chaotic carrier the more diflicult will be to decode the message without the appropriate 
receiver because of the high frequencies and the large number of degrees of freedom of the chaotic carrier. It is 
worth noting that the chaotic carriers in semiconductor diode lasers provide a broad spectrum in which the message 
can be hidden. With the advent of high speed telecommunications and advancements in computer processing [14] 
technology, the demand for faster communication is fulfilled with these laser systems which are increasingly used in 
major applications. 

Researchers are now looking at ways to exploit lasers with chaotically fiuctuating signals, to add an extra layer of 
privacy to messages sent over fiber-optic lines. Such optical systems will generate a chaotic carrier with considerable 
high dimensionality which will enhance the security level together with very high transmission rates. These higher 
dimensional nonlinear devices exhibit nonlinear dynamical behavior with fast irregular pulsations of the optical power, 
or wavelength hopping, with bandwidth ranging from a few giga- hertz to tens of giga- hertz, large correlation dimension 
of chaotic carriers [15,16] preventing linear filtering and frequency- domain analysis, complexity and unpredictability 
of the chaotic carriers, compact and integrated devices, high dimension due to the delay-induced dynamics [17] etc. 
In recent years, experiments were conducted with optical chaos based communication involving fiber lasers [18] and 
semiconductor lasers [19]. Then Argyris et al [20] in 2005 established high speed optical communication of IGbps 
over a commercial fiber-optic channel. This was a pivotal point in improved communication over other commercial 
and common used modes of data transmissions. The chaotic optical setup generated an optical carrier to encode and 
transmit a message over a communication distance of 120 km achieving a robust transmission system, immune to 
perturbations and channel noise. A major requisite for optical cryptography is to convert the message to be sent into 
an optical signal. 

In order to realize secure communications using semiconductor lasers, the emitter must exhibit a chaotic regime 
which is subjected to chaotic oscillations by either injection from another source [21] or by refiection from an external 
mirror [22]. This can be achieved by when the emitter is either subjected to an optical or an electro-optical feedback. 
Our system is used as an optical feedback, illustrated in Fig.l, with emitter and receiver semiconductor lasers 
unidirectionally coupled. The setup consists of semiconductor lasers which are biased above a threshold level. These 
are then subjected to optical feedback using external mirrors. The receiver operates as the emitter under similar 
conditions in a closed loop or under CW when decoupled by open loop scheme. A piezo transducer is employed when 
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FIG. 1: (color online) Schematic diagram depicting communication scheme with two semiconductor lasers subjected to optical 
feedback with the emitter and receiver lasers being coupled unidirectionally. 

working in a closed loop scheme, so as to match the external cavities. The optical signal thus produced is fed into a 
laser that passes it along the laser beam. Then chaotic synchronization occurs between the two spatially separated 
semiconductor lasers wherein the receiver side tunes its oscillations to that of the emitter but suppressing the hidden 
information. At this stage, the message is encrypted into the chaotic carrier. 

Synchronization basically implies that the receiver can reproduce the exact chaotic irregular orbit of the emitter, 
i.e. the transmitter and the receiver both generate the same chaotic signal by reproducing the irregular time evolution 
of the same chaotic carrier, which traces a chaotic orbit at the emitter end. Then the chaotic carrier wave within 
which the message is embedded is transmitted from the emitter laser to an almost identical laser, which forms the 
receiver section. The output of the receiving laser is tuned to match the incoming signal from the transmitter. During 
decryption, the laser in the receiver section only recognizes the chaotic mask of the incoming information-bearing 
signal and extracts the originally transmitted message. It extracts the message by subtracting the emitter laser beam 
(the input), which contains the message from the identical carrier (output) generated at the receiver. 

Considerable research [23-25] shows how efficient optical feedback laser systems are. In 1997, the authors [26] 
applied synchronization of the optical chaos from external reflection technique based on optical feedback for a robust 
cryptographic application by chaotic shift-keying, which is resistant to external perturbations and noise. There have 
been other signiflcant successful attempts in chaos based optical communication and its application in cryptography 
[27-34]. However, there are other issues to be addressed regarding the reliability and efficiency of optical cryptography, 
concerning parameter mismatch, noise, dispersion and modulation index of chaotic carriers. Note, that the parameter 
mismatch has to be prevented since the correct decoding at the receiver section intrinsically depends on the tuning 
accuracy of the master and slave semiconductor lasers [35]. To facilitate decryption the chaotic system possesses 
identical parameters and operating conditions. If there is any mismatch, the devices will not be synchronized and 
hence will not fulfill the goal of cryptography. Moreover, the scheme has to overcome the hurdles imposed by the 
high-frequency together with the high dimensionality of the chaotic signal. 

Another important issue is the noise produced in the lasers or in the channel communication. Such noise will 
hamper the accuracy of the synchronization due to which the synchronization succumbs to destabilization. Then, the 
dispersion factor of a chaotic laser is also an area of concern. The chaotic signal is highly sensitive to dispersion, since 
it has a broad line width in the order of Giga Hz in comparison to the output radiation of a solitary laser which has a 
line width of only few Mega-Hz. Therefore, due to the increased line width, the fiber communication links should be 
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shorter than the traditional ones [36] to facihtate propagation by chaotic semiconductor lasers. Lastly, the modulation 
index of the encryption should be small to achieve enhanced security. 

In this paper, we have endeavored to use synchronized optical feedback induced chaotic dynamics for color image 
cryptography. The latter uses an external cavity semiconductor laser on the transmitter side and a semiconductor 
laser with optical injection on the receiver side. Image cryptography is achieved by forming different masks derived 
from preprocessed chaotic sequences arranged in the form of an image whose elements are chaotic sequences. 

In the case of the two lasers, one should note that both systems are highly nonlinear and infinite dimensional in 
nature, even if one may only simulate a finite dimensional approximation, in a computer. The respective dynamics 
are also very sensitive to initial conditions, but it is possible to synchronize the two systems with a proper optical 
feedback mechanism. As we know, the synchronization of two continuous spatiotemporal PDE systems are quite a 
new subject of investigation. So, in our study, we have found interesting results which may motivate a new line of 
research in cryptography. 

It is remarkable that, although both systems are infinite dimensional, they may easily synchronize with different 
initial and boundary conditions, using proper one way couplings. The experimental setup can be implemented as 
illustrated in the paper (Fig. 1) and the synchronization phenomenon between two lasers can then be used for practical 
purposes. The cryptographic scheme is just an example of the possible applications in the field of communication, 
which profits from the features of the system. The differences between any synchronized chaos based cryptosystem 
and our scheme are the following [1]: 

1. The driving and the response systems are highly nonlinear and infinite dimensional in nature. Once the two 
systems synchronize, they remain synchronized forever. If we have synchronization at t = tmax = 5 (dimen- 
sionless unit), they will remain the same for any ti > tmax- Therefore, both the sender and receiver can start 
collecting data from any time after tmax- They can choose the data in an infinitely long time interval. 

2. We have three variables in our main PDE system [Eq.(l)] viz. the forward field the backward field 
E_i and the normalized carrier density function TV. The keys used to encrypt /decrypt the image file are the 
variables E±i^E±2 and Ni,N2 respectively. The interesting thing is that the driving signal from sender to 
receiver constitutes a one way coupling and contains only the electric fields (eq 4) of the driving system; it does 
not contain any information about the carrier density. In this equation, Ni as well as N2 are not transmitted. 
Hence, the set of keys are 

Private/ secret key for transmitter Ni 
Private/ secret key for receiver N2 
Private key for both parties tmax 

Public key for the transmitter E_^i,E_i 
Public key for the receiver Ej^2^E-2 

3. Therefore, any third party attempting to hack the signal sent from the transmitter to the receiver, will only 
get information about the electric field E but not about the carrier density, which is used as our secret keys for 
encryption. They can try to reconstruct the phase space on the basis of the data obtained from the time series 
of the electric field E, but it is impossible to predict the nature of such a infinite dimensional PDE system. 

These considerations motivated to exploit the potential of such a complex system in secure communication. The 
pre-requisites for any cryptographic application are met by our scheme: the system is high dimensional, guaranteeing 
an enormous key space, and is highly sensitive to minute perturbations, since it is spatio-temporally chaotic. This 
justifies consideration of our method for cryptographic applications. The system which is used to devise the algorithm 
for cryptographic application generates an infinite key space. 

• So after synchronization, and after the time tmax^ a state is chosen, from which two sets consisting of infinitely 
many numbers, and F^, are generated. 

• A key should be large enough that a brute force attack (possible against any encryption algorithm) is infeasible, 
i.e. would take too long to execute. According to Shannon's work on information theory, in order to achieve 
absolute secrecy, it is necessary for the key length to be at least as large as the message to be transmitted and 



[1] Detailed definitions of the quantities which we briefly mention here are given in the following sections. 
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should be used once. This algorithm is called the One-time pad, also known as the Vignere cipher. Since longer 
symmetric keys require exponentially more work under brute force search, a sufficiently long symmetric key will 
prevent this type of attack. In our scheme, the key length is as large as the message which is transmitted [Eqs. 

(6) , (7)]. Through the modulo operation, the length is restricted to the range [0, M] or [1, M+1]. These processed 
keys are of the same length as that of the message and are used exactly once and therefore this organization is 
theoretically unbroken [4]. So, the advantage of the Vignere cipher is reaped by our scheme. Then in order to 
encode an image of length M x A/", we convert the set and Y'^ into a chaotic sequence Ki and K2 respectively 
containing the M x N elements as that of the plain image P of dimension M x N. Each element of the set 
KiandK2 are integers in the range [0, M] or [1, M + 1]. This conversion is through eq. 6 and eq 7. For example 
as illustrated later, with a 512 x 512 sized colored image P we have restricted the key space. We obtain a set 
for P by transforming into an array of integers in the range [0,512] from the processing steps in Eqs. (6) and 

(7) . 

• Although, chaos is an irregular motion, it is deterministic and therefore the original image can be completely 
recovered if the secret keys are exactly known. But in our scheme, two of the four secret keys, the carrier density 
and tmax are not transmitted over the network and only the algorithm is assumed to be public. As a result, our 
scheme adheres to the Kerchkoff 's principle for cryptographic algorithms. 

The rest of the paper is organized as follows. In Section. II, we define a set of spatiotemporal partial differential 
equations (PDEs) that describes the dynamics of semiconductor lasers. The system exhibits chaos for some particular 
values of the parameters. In Section. Ill, we have implemented the synchronization phenomena in a set of coupled 
PDEs with linear feedback coupling. The physical significance and the results are described in detail. Section IV 
presents the encryption scheme based on the synchronized system of PDEs. We study the security analysis of the 
proposed cryptosystem with some statistical analysis in Section. V. Finally, Section. VI reports our conclusion. 



II. SPATIOTEMPORAL SEMICONDUCTOR LASERS AND THEIR CHAOTIC PROPERTIES 

A spatiotemporal chaotic system is a spatially extended system that exhibits spatiotemporal chaos (STC), i.e. 
higher disorder in both space and time. The disorder in space can be related by rapid decay of spatial correlations 
whereas the temporal disorder can be measured by the positive Lyapunov exponents. Examples include the one 
dimensional coupled map lattices CML, which are discrete-time and discrete- space dynamical systems consisting of 
nonlinear maps acting on the lattice sites, coupled with each other. Because of the nonlinear dynamics of each local 
map and the diffusion due to the spatial coupling, a CML exhibits STC [37]. Another example of a system exhibiting 
spatiotemporal chaotic regime is the semiconductor laser system. The advancements in telecommunications have 
opened up a new arena in the study of nonlinear systems which include the high nonlinearity demonstrated by laser 
systems. Individual chaotic systems are often marred by the problem of rapid degradation of the chaotic dynamics 
in finite computing precision [38]. It is in this area that spatiotemporal systems have found an edge, for such chaotic 
systems have a sustained periodicity [36] with a performance superiority in cryptography [39] . 

Investigations with optical systems have revealed that the dynamics of laser systems with optical feedback depends 
on the size of the external cavity. The semiconductor lasers have been found to oscillate with several longitudinal 
modes simultaneously [40]. Such multimodal behavior increases the dimensionality of the rate equations and has 
been extensively used in applications, such as fiber couplers or compact discs. This is because semiconductor lasers 
have short external cavities. To realize this, one places an external refiector which feeds back a fraction of its delayed 
electromagnetic output into the active layer of the semiconductor laser. If the distance between the refiector and 
the output facet of the semiconductor laser is L^xti then the round-trip time in the external cavity is 2Lext/c. A 
convenient model of this configuration should be imperatively made in the field representation, because the phase of 
the coherent laser radiation is strongly affected by the delayed feedback. Therefore, the system we consider takes the 
following form [41]. 

±^ + 1^ = ' (r ) - + - [(iv - 1) - »ivi E,. 

where E± are the forward (backward) electric fields normalized by the electric pump Eq ~ 10^^ v/m and N is the 
carrier density normalized by A^o- The space and the time variables are normalized respectively by Lq ~ 1/i m 
and T ~ 0.01 ns. The system of Eqs. ([T]) describes the weakly nonlinear dynamics of slowly varying two counter 
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propagating longitudinally high-frequency (hf) optical fields (traveling waves) interacting with low- frequency (If) 
density perturbations associated with the electron-hole (total) charge carriers in the active semiconductor lasers. In 
the case of a conservative system, the formation of optical envelope solitons through the nonlinear interactions of 
such waves is most important in the context of turbulence (in which energy is transferred to few interacting wave 
modes) as well as in stable optical pulse propagation. Here the amplitude of the optical wave envelope is modulated 
by the electrostatic small but finite amplitude electron-hole density fluctuations in which ni represents the refractive 
index of the active layer, Dp the diffraction which provides higher-order dispersion. All other parameters and the 
corresponding numerical values for exhibiting chaos are given in Table I. Physically, since the coupling of the two laser 
strips, which is provided by the overlapping evanescent optical fields and by diffusion of electron-holes, can be changed 
by varying the inter-element distance, the dynamical behaviors of the twin- strip lasers can show a transition from 
order to spatiotemporal chaos of the optical field intensity [41]. Moreover, since the optical properties are strongly 
affected by the local optical fields and the properties of the charge carriers along the longitudinal directions, the 
dynamical behaviors of them are to be determined in a self-consistent manner as described above by the system of 
three PDEs. We numerically solve the Eqs. ([T]) by Runge-Kutta method with initial conditions as [42,43] 



TABLE I: Laser Parameters for computations 



Parameter Value#l 

L (Cavity length) 275 fim 

w (Strip Width) 5 /xm 

d (Active layer Thickness ) 0.2 /xm 

Ri (Reflexivity of the front mirror ) 0.15 

R2 (Reflexivity of the back mirror ) 0.95 

R3 (Reflexivity of the external cavity) 0.95 

ni (reflexive Index of the active layer ) 3.5 

,trim=O.Oin Oin Oin -1.5in a (Linewidth enhancement gain) 5 

a (gain coefficient) 1.5x 10~"^^cm^ 

No (transparency concentrations) 0.67x 10"^^cm^ 

A (Laser wave length) 8000 nm 

T (nonradiative recombination time) 5ns 

Text (Round trip time for the external cavity) 830 ps 

Lext (External cavity length) 12 cm 

Df (diffusion coefficient) 30 cm^/s 

Dp (Refraction Coefficient) 10~^ m 

r] (Injected efficiency) 0.5 

J (Injected current) 80/95/110/130 mA 

aint (Internal loss) 

s (Surface recombination velocity) lO^m/s 

13 (output scaling factor) 1 

r (Confinement factor) 0.5 

ric (refractive index of cladding) 3 

rieff (Effective index) 3.5 

A^o (Career density at transparency) lO^^m""^ 



E±{z,0) = E±o^E±icos{kz)/P, N = -N cos{kz)/ (2) 

where k is the wave number of modulation, /3 is a suitable constant taken as 1/500 to ensure that the perturbation 
is sufficiently small, E±o = 2 and £^±1, N are constants of the order of unity. The numerical results are shown in Fig. 
2 after the end of the simulation at t = 100. We increase the injected current from J = 80 to 130 mA with external 
cavity length L^^t = 15 cm and round-trip time r = 830 ps. We choose the injection efficiency 77 = 0.8 to assume that 
80% of the carrier reaches the active region. From Fig. 2, we find that the forward as well as the backward fields 
are chaotic in nature along with the density perturbation, and we can optimize the maximum electric fields to obtain 
values greater than \E± \ = 20 for both waves. 
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FIG. 2: (color online) Numerical solution of Eqs. ([T]) for forward and backward electric fields (see the upper two panels) as 
well as the carrier density (lower panel) exhibiting spatiotemporal chaos. 

III. SYNCHRONIZATION OF SPATIOTEMPORAL SEMICONDUCTOR LASERS 

In this section we will synchronize the two nearly identical semiconductor laser systems [c.f. Eqs. (IT])] which exhibit 
chaos as shown above. To this end, we use linear coupling terms with a coupling coefficient e. The former behaves 
like a damping rate of the optical fields in order to balance or manifest the nonlinear laser loss or gain, and can be 
controlled through the variation of the coupling strength e. Thus, we have the following system of equations where 
the suffix 1 is used to denote the driver, and 2 the response system. 



± 



dz 







d^E±, 


dt 




dz^ 


dNi 




\ d'^Ni 


dt 


-m: 


I dz^ 



iWLeE±i + k[(N-1)- iaN] E±i, 

iVi - aTE^o {Ni - 1) (l^+il' + |^_i|') , 



AT 

iVo 



T 

T 



(3) 



and 



dE 



±2 



dz 



dE±2 




d^E±2 


dt 




dz^ 


dN2 




\ d^N2 


dt 




1 dz^ 



- iWLeE±2 + k[{N-1)- iaN] E±2 + e {E±i 
AT T . 



E±2) , 



-No 



T 



aTEl{N2-l) (\E+2\^ + \E_2\^'^ 



(4) 



where n = VLcN^a^e = (Lc^Lc = Tc/ni. Thus, the two coupled multi-mode chaotic semiconductor lasers given by 
Eqs. ^ and Q describe the driver and response system respectively to demonstrate chaotic communications in 
high-dimensional systems through synchronization. This will enable the transmission of multiple messages through 
a single or multiple scalar complex field. The synchronization of the two systems ([3| and Q, is achieved through 
the inclusion of coupling terms ex e which, on the other way, describe non-radiative laser loss or gain by the driver 
or response field due to power reflectivity. It is, therefore, of importance to see whether synchronization of such 
high-dimensional chaotic systems is robust to noise and/ or error signals for transmitted messages. Using the similar 
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FIG. 3: (color online) Synchronization errors of coupled spatiotemporal semiconductor laser systems. 

initial conditions along with E±^N and their derivatives vanish at the boundaries, we integrate the coupled systems 
(|3| and Q numerically. The simulation results are displayed in Fig. 3. These basically represent the synchronization 
error between the driving laser ([3| and the response one Q respectively at the most chaotic state. The initial values 
for numerical simulation are given in Eq. (2) with ^+io = Ni = 0.2, ^+ii = 1, E-io = 3, E-u = 1. We find that the 
corresponding errors for the forward and backward electric fields are of the order 10~^^ and for the carrier density the 
difference is of the order 10~^^. So, at the synchronization state the system ^ can be considered as the transmitter 
and the system Q as the receiver one. 

IV. OPTICAL CHAOS AND APPLICATION TO IMAGE ENCRYPTION 

This section focuses on the application of the synchronized space time laser systems to digital cryptography. Image 
encryption basically denotes the realignment of the pixels to different positions with a change in the statistical value 
of the pixels. This permutation yields a disordered version of the image which is totally incoherent. In order to 
realize image cryptography, it is imperative to devise a transmitter and receiver section. It is in this arena that 
synchronization is utilized. For cryptographic encoding, the transmitter for our scheme uses system (3) and the 
receiver is formed from system (4). Both of them are certain to chose values from the time series E±i, Ni and E±2^ 
N2 respectively after synchronized space and time for the proposed cryptographic scheme. The proposed method is 
simulated in Matlab 7.0 environment. 



A. Generation of chaos sequence 

A colored digital image P(M, D) is basically composed of three sets of two-dimensional (2D) matrices containing 
integers in the range of — 255. Each matrix contains M number of rows and TV number of columns and has three- 
color planes, D namely Red (R), Green (G) and Blue (B) which are concatenated along the third dimension D. Each 
picture element of the three-dimensional (3D) matrix in M x A/" x 3 format is actually a triplet known as a voxel. It 
is analogous to its 2D monochrome image, the pixel. In this work, we would henceforth refer a picture element for a 
colored image as a voxel. Due to the availability of fast processing and high bandwidth technology, it is feasible to 



8 




FIG. 4: (color online) Illustration of the encryption scheme: (a) Encryption by transmitter- T . Symbols :- A - Chaotic data set 
y^, B - Transformation along the color planes to form CSI^ , C - mask CSI^ ^ D- Shuffling operation in encryption, E - The 
shuffled image T, F - Chaotic data set I2, G - Transformation along the color planes to form CSI^, H - the second encryption 
mask CSI^ ^ I - Transpose operation, J - The transposed mask CSI^ , K - differential mask CSI^ . 

employ encryption directly onto the 3D image instead of pre processing it to an indexed 2D monochrome image. This 
reduces considerable computational costs. Hence, the plain image P containing len = 3MN voxels expressed by 

P = {pi =Pi,P2 =pf,P3 =pf,P4 =P2^P5 =P2^P6 = , . . . , Pi , . . . , PZen =Pmn}- (5) 

B. The algorithm for image encryption 

In digital image cryptography, Confusion and diffusion are considered as two distinct stages, both requiring image- 
scanning to obtain pixel/voxel values. In our method, the original message undergoes a complete randomized permu- 
tation of the positions in the diffusion phase and a drastic statistical change in the confusion phase. The proposed 
method eliminates the necessity of duplicated scanning thereby enhancing the encryption speed. At first, the image 
is partitioned into blocks of voxels. Then STC is employed to shuffle the blocks with a simultaneous change in the 
voxel values. Also, an efficient method for generating pseudorandom numbers from STC is suggested, which further 
increases the encryption speed. Theoretical analysis and computer simulations confirm that the new algorithm is 
highly secure and is very fast for practical image encryption. 

The stream cipher mode of encryption has been implemented at the bit level with RGB image. We employ the 
advantages of ergodicity and sensitivity to initial conditions of synchronized chaotic system which leads to the 
confusion and diffusion properties required for secure cryptosystem. The entire cryptography scheme is illustrated in 
Figs. 4 and 5 and explained as below. 

Preprocessing Stages' - 

1. :- The system under consideration is subjected to a maximum number of iterations tjnax- An initial synchronized 
state is selected to extract the chaotic sequences from variables E±i^ Ni at the transmitter end and £^±2, 
N2 at the receiver section. Then, the i -th chaotic data for permutation of the elements of the RGB image P 
of size 3 X M X N is obtained from a synchronized data set = {yl^yl, . . . , ?/^} and = {yl^ • • • 5 
respectively from £^±1, Ni and are preprocessed since the elements are real numbers. The scheme employed 
is a duplex encryption technique. Hence, the two different data sets of keys, namely = {kl^k^^ . . . ^k^^} 
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FIG. 5: (color online) Illustration of the decryption scheme: (b) Decryption in Receiver section - R. Symbols :- A - Chaotic 
data set y^, B - Transformation along the color planes to form CSI^ , C - mask CSI^ ^ D- Ant i- shuffling in decryption, E - The 
shuffled image T, F - Chaotic data set I2, G - Transformation along the color planes to form CSI^, H - the second encryption 
mask CSI^ ^ I - Transpose operation, J - The transposed mask CSI^ , K - differential mask CSI^ . 

and = {/c^, /c^, . . . , k'^} that are used to encrypt the image are generated through the simple algebraic relation 

kl ^ integer(abs(10^ x \yl\))Mod{M + 1), (6) 

and 

k'^ ^ integer(abs(10^ x \y'^\))Mod{M + 1). (7) 

Hence, i^^and K^e [0, M] 
2. :- 

• Next we select a subset of integers from the processed data set and i^^from Eqs. (6) and (7) into 

respectively which will now be an array of M x data elements with M rows and N columns and 
further process these through the following rules 

K^' ^ reshape(i^^ , M, A/") , (8) 

K^' ^ reshape(i^^ M, N). (9) 

where reshape means to change the size of and to M x N array, from which the processed elements 
are extracted in a column wise fashion. The result of this operation will be an array containing an M x 
number of elements. 
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a 



b 



FIG. 6: (color online) The masks used in two phase cryptography: (a) Mask generated from first chaotic sequence (b) Mask 
generated from second chaotic sequences. 

• and can be represented as an image. In our scheme, we have tiled the array which have resulted 
from Eqs (8) and Eqs (9) along the three color planes by using the command repmat(X,[P Q]) in Eqs. 
(10) and (11). This will produce a large array consisting of P x Q block of array where X may be a 
multi dimensional array. Since, in this operation we generate an 1 x 3 tilings of copies of and 
respectively, we denote it as equivalent to an image matrix which we label it as the Chaotic Sequence 
Image CSll^ j^j^, CSlj^ j^ jj respectively. 

CSll^^N^D ^ repmat(i^i', [1 1 3]), (10) 

CSlj^^N^D ^ repm^t{K^\ [1 1 3]). (11) 

Here, CSI^ will serve as the permutation mask [Fig. 6 (a)] in the diffusion stage and CSP [Fig. 6(b)] is 
used during the next phase of encryption. Each element of these phase masks are composed of chaotic 
sequences and therefore are completely randomized themselves. 

3. :- The CSP is transposed into CSP' as 

^^4%,,, (12) 

where by r = {1, 2, 3, ... , A^}, c ={1,2,3,..., M} and d = {1, 2, 3} Then the statistical difference between 
CSI^ and CSP is computed as 

CSI^ ^ abs(CS'/^ - CSI'^\ mod (M + 1)). (13) 

Encryption:- 

Position Shuffling :- Diffusion phase 
Step A:- 
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Sort the elements of CSI^^^ in ascending order in accordance to their voxel values as a vector where 
r = {1,2,3, ...,M}, c = {1,2,3, ...,7V} and d = {1,2,3}. Let this be known as CSll[^^^. Each of the ele- 
ments of CSI^ form the indices for position scrambling of P. 

Step B:- 

Shuffle all the M columns and 7V rows of the original image P{p^,p^ ,p^) with CSI^ into the intermediate 
transformation matrix Tr^c,d by the following simple process 

Permutation along the rows- Here we extract the pixel value of along its rows indicated by the symbol: to form the 
keys for row shuffling 

Tr,c,d Pkj,c,d- 

Next, the columns of the intermediate resultant image is also replaced by CSI^ . Permutation along the column is 
shown by a similar procedure where the pixel values along the column are extracted to form the keys as under 

Tr,c,d ^ Pr,kl^^,d- 

Statistical Alternation:- Confusion phase 
Step Cl- 
in this additional phase, the matrix obtained in Eq. (13) is used to change the binary information of the im- 
age by the XOR operation represented by the symbol 0. Thus, the following step is executed to get the cipher 
image CI 

CI ^T^CSI^. (14) 

The operations discussed here occur along the three-color planes producing the final encrypted image CI. This is 
the cipher image which is transmitted along with the initial conditions as the secret key. 

Decryption:- 

Step D:- 

The same data set and Y'^ at the receiver end is preprocessed as in preprocessing steps 1 to 3. Then we 
perform anti-XOR and reverse permutation operation as in step A to C to recover the original image P. These 
procedures are lossless and reversible in nature and thus, without any overhead, we can correctly decipher the image. 



V. SECURITY ANALYSIS AND SIMULATION RESULTS 

This section uses the algorithm for encrypting a 512 x 512 colored image and examines its security. The bands 
in Fig. 7 (c) results from the abrupt changes in the diffusion process due to the inherent dynamics of the nonlinear 
system under consideration. It is noteworthy, that a row or a column may be permuted more than once which 
may be allowed since the computational cost in linear check of a duplicate value in the data set is higher than the 
ultimate effect of eliminating redundant position values. Figure 7(d) depicts the complete randomized orientation of 
the voxels making this phase totally incoherent from the original image. 
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FIG. 7: (color online) The encryption process: (a) The original color image (b)The mask for permutation of original image (c) 
The shuffled intermediate image (d) The final cipher image. 



A. Effectiveness of the scheme 

The proposed method correctly recovers the original message. This is shown by a metric known as Mean Square 
Error (MSE) given by 



A correct recovery of the image will yield MSE = [MSE^ MSE^ MSE^] = [0 0] for all the three-color channels 
separately. It is observed that the recovered image highly differs from the original with a minor change in the key set 
as explained subsequently. 



B. Key space attack and sensitivity analysis 

For the numerical simulation of the system (3) and (4), we choose the initial conditions as 




M N 



(15) 



r—1 c=l 



E+i{z,0) 



E+io + E+ii cos{kz)/L, E+2{z, 0) 



£■+20 + £'+22 cos{kz)/L, 



(16) 



E-i{z,0) 



E_io + £-11 cos{kz)/L, E_2{z, 0) 



£_2o + £-22 cos{kz)/L, 



(17) 



Ni = -Ni cos{kz)/L, N2 = -N2 cos{kz)/L. 



(18) 



where 
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^+10 = Ni = 0.2, E+n = l;E_io = 0.3, ^_ii = 1.1; ^+20 = 0.5, ^+22 = 1-2; ^-20 = 0.6, ^_22 = 1-3, iV2 = 
0.5,1/ = 1/500, e = 2.5 (in arbitrary units). 

The total number of different keys that can be utihzed in a process is termed as the key space. The 128 
bit paradigm is a widely used concept in the AES cryptography scheme which is commercially used nowadays. It 
is a block cipher mode under the symmetric cryptography scheme. This falls under the category of traditional 
cryptographic technique including DES, RSA, DSA, hash functions, one time pad etc. However, our method based 
on a chaotic symmetric cipher stream mode of encoding the message, is an offshoot of the traditional methods of 
cryptography. In order to render high security, it is imperative that the key space be vast enough to make any 
brute force attack ineffective. Also, to ensure a large divergence of a chaotic trajectory from the initial condition, 
the number of iterations should be relatively large. Both these issues have been addressed. The key space with its 
keys are generated after a considerable number of iterations which consists of the private and the public variables 
viz. E±i^E±2 and TVi, 7V2. So, due to these keys the system has a pretty exhaustive key space which is large enough 
and as such makes our system resistant to brute force attacks. Each key has a value which lies within — 512 for 
a colored image of size 512 x 512. These are used to encode the color values of the 9 bit RGB image after tiling 
operation along the 3 color planes according to Eqs. (8) and (9). 

An encryption technique is good and efficient if the cipher image possesses sensitivity to the changes in the cipher 
keys, together with a combination of sufficient large key space, in order to render the system immune to major brute 
force attacks. 

A wrong data set in any one of the stages will result an incorrect decrypted image as shown in Figs. 8 and 9. The 
sensitivity to keys is tested on an RGB image of size 512 x 512 x 3 through the following test cases :- 

1. Testl:- Let us examine the effect of alternation of one of the secret keys at the receiver section. The original set 
of keys derived from variables used by the transmitter are £^±1 and A^i from Eq. (2). E±i =0.2 is one of the 
initial condition of the system which is also of the the public key. However, the initial condition values are 
never transmitted over the network. Now, let us assume that an attacker intercepts the communication network 
and self-imposes itself as the transmitter. Since, the initial conditions are never transmitted over the network, 
the invader attempts to make a guess about the initial conditions. In this process, let the guessed value be 
assumed to be E±i = 0.200001. 

But, the key set used for decoding at the receiver section will not match with the ones used in encoding the 
message as the attacker has its own wrongly generated keys. So, the effect of this tampering is reflected in the 
incorrect recovery of the message at the receiver section in Fig. 8 (f). 

2. Test2:- The keys in the permutation matrix Eq. (6) and diffusion matrix Eq. (7) are altered through the 
following formulae 



The key space increases in dimension by further increasing the dimensionality during the processing stage in 
eq 6. In this part, we show that the keys so obtained are highly sensitive to slight change. Hence, if somehow 
or in a remotely possible case, another permutation matrix is generated by a third party who has gained an 
unauthorized access, the result will be an incorrect recovery of the message. So, the position permutation matrix 
is changed in Eq. (19) to display this effect of tampering without any assumption of changing or reducing the 
power of 10. Thus, the first key changes from 190 to 33, second key changes from 294 to 322, third from 398 to 
362 and so on. The outcome is demonstrated in Fig 9 which also has high MSE values. Such high MSE values 
show that the image obtained by the receiver is not even close to the original and is highly uncorrelated from 
the original from which no visual interpretation can be deduced. Therefore, the receiver will come to know that 
the key set has been tampered. By Eq. (20) we have shown that the keys used in diffusion phase are equally 
sensitive to minute disturbances on being tampered by an attacker. So, without any specific reason for the 
choice of the power of 2 which is raised over 10, we have modified the key space with a reduced dimension of 
10 in eq 20. The outcome is shown through the MSE values listed in Table II. Needless to say, through these 
exhaustive statistical tests, the sensitivity of the key sets to perturbations is validated. The effect of these are 
demonstrated in Fig. 9f where the decrypted image is totally incoherent from the original and there is no way to 



kj ^ integer(abs(10^ x \yl\))Mod{M + 1), 



(19) 



k'^ ^ integer(abs(10^ x |^^^|))Mod(M + 1). 



(20) 
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FIG. 8: (color online) Test result of statistical Test 1: (a) Chaotic phase mask CSI^ for the permutation phase with the 
correct keys and initial conditions (b) Second chaotic phase mask CSI^ obtained from the correct keys and initial conditions 
(c) Correct decrypted image (d) First chaotic phase mask for the permutation stage obtained from the incorrect keys and 
tampered initial, trim=O.Oin Oin Oin -1.5in conditions (e) Second chaotic phase mask obtained from the incorrect keys and 
tampered initial conditions (f) Incorrectly decrypted image 

predict any sort of information about the original message from the incorrect recovered result. The preprocessing 
of the key sets along with the transpose operation of the second phase mask which is then used to to form the 
mask image for the confusion phase adds enough stochastic properties to render the cipher image totally obscure 
from the original image. 



TABLE II: Mean-square error for Test 1 and Test 2 



MSB 


R G 


B 




Testl 
Test2 


l.Oe+003 X 9.0296 1.0e+003x c 
l.Oe+003 X9.1335 l.Oe+003 X( 


B.6193 1.0e+003x8 
5.7111 l.Oe+003 x 


.6450 
8.5998 



The MSE value listed in Table II shows that incorrect decryption results in a very high MSE between the original 
and decrypted image. The values are not even closer to zero since the recovered output is remotely unrelated to the 
transmitted original image. Moreover, from Figs. 8 and 9 we can infer that it is almost impossible to reconstruct 
the original image when the secret keys and the data sets are tampered. This further reinstates that the process is 
immensely sensitive to small variations to the initial conditions. 

C. Differential attack resistance 

A good cipher should be immune to the influence of a change in a pixel. The rate of change of the number of 
pixels of the cipher image CI when the plain image P is altered by just one pixel is quantified by a metric known as 
NPCR. It is defined by 



NPCR = '^Yr^ J' X 100, (21) 



15 



FIG. 9: (color online) Test result of statistical Test 2: (a) Chaotic phase mask CSI^ for the permutation phase with the correct 
keys (b) Second chaotic phase mask CSI^ obtained from the correct keys (c) Correct decrypted image (d) First chaotic phase 
mask for the permutation stage obtained from the incorrect keys (e) Second chaotic phase mask obtained from the incorrect 
keys (f) Incorrectly decrypted image. 

where D(r^c) represents the change in the picture element due to encryption for a monochrome image. 

D{r^ c) = when P{r, c) = C/(r, c)else 

D{r^ c) = 1 when P{r^ c) ^ C(r, c) 

In our example, the change in a few voxel values before and after encryption are listed below to demonstrate the 
dispersion of voxels through the scheme 

P(l, 1, 3) = 153, C/(l, 1, 3) = 218, 

P(2, 2, 3) = 148, C/(2, 2, 3) = 74, 
P(50,50,3) = 152,C/(50,50,3) = 113, 
P(100, 100, 3) = 79, C/(100, 100, 3) = 238, 
P(256, 256, 3) = 42, C/(256, 256, 3) = 123. 

This shows that the resulting cipher image has its picture elements totally changed. Next quantifiable measure of 
diffusion properties is the Unified Average Changing Intensity (UACI) which determines the average intensity of the 
differences between the pixel values of the original and encrypted image. It is expressed by 



M X TV ^ 255 ^ ^ 

r=l,c=l 

The results of NPCR and UACI in percentage is shown in Table III. From it we can deduce that the scheme has a 
high value of NPCR along with a satisfactory value of UACI. Thus, higher the value of NPCR, more is the system 
resistant to differential attacks. This validates the efficiency of the proposed scheme. 
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TABLE III: Result of differential attack: NPCR and UACI performance metric 





RGB 


NPCR 
UACI 


99.61 99.62 99.60 
14.47 6.97 10.35 
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FIG. 10: (color online) Histogram: (a)-(c) original image, (d)-(e) cipher image. 



D. Histogram analysis 

The distribution of the pixels after being scrambled or manipulated can be studied with the help of a chart which 
represents the distribution of the pixels in the range — 255. For a gray image having gray level Vq it is represented 
by a function hist(/g) = n^, where Iq denotes the q^^ gray level and Uq are the number of pixels in an image. In the 
diffusion phase, the positions of the image elements undergo shuffling. This does not alter the statistical information 
of the original image. However, the additional layer of security will disguise the desired information, as reflected in the 
histogram of Fig. 10. The intensity of the encrypted image greatly differs from the original image, thus an attacker 
will be unable to infer any statistical information required to decode from the scheme. From the result of intensity 
distribution, it can be ascertained that the scheme also possesses good confusion properties. 



E. Correlation analysis of pairs of adjacent pixels 

Statistical analysis of the cipher image can be performed by another excellent metric, the correlation coefficient. 
This shows whether there is any occurrence of an association between the guessed secret key and the observed 
information of the cryptosystem and forms the core information for cryptanalysis. The covariance cov(x, y) between 
a pair of pixel values present in the subset x and y in sl gray image is defined by 



cov{x, y) = E{x - E{x)){y - E{y)), 
Then, the correlation coefficient Vxy is given by 



_ Cov{x,y) 
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FIG. 11: (color online) Scatter plot of correlation coefficients of adjacent pairs of pixels along the horizontal direction: (a) 
Correlation coefficient of the original image in the Red color plane (b) Correlation coefficient of the original image in the green 
color plane (c) Correlation coefficient of the original image in the blue color plane. 



where E{x),E{y) denotes the mean; D{x)^D{y) stands for the variance between the pixels. Correlation coefficient 
values of a random sample of 1000 voxels, placed adjacently along the horizontal, vertical and diagonal directions, 
have been computed for each of the three-color channels and enumerated in Table IV. From the values, it can be 
observed that the correlation of cipher image is almost zero. For an efficient encryption methodology, it is imperative 
that the correlation values of adjacent pixels be minimal for the encoded image. Our scheme fulfills this criteria. Fig. 



TABLE IV: Correlation Coefficient for each color band 





Original 


Cipher 


R 


G 


B 


R 


G 


B 


Horizontal 

Vertical 

Diagonal 


0.9714 
0.9592 
0.9621 


0.9756 
0.9661 
0.9690 


0.9730 
0.9610 
0.9640 


-0.0002 

0.0429 

0.0305 


-0.120 
0.0539 
0.0270 


0.0207 
0.0793 
0.0281 



11 displays the spatial distribution of the randomly selected voxels from the three-color planes placed adjacently along 
the horizontal direction. The graphical result emphasizes that there is hardly any correlation between the original 
and the distorted version of the image. 



VI. CONCLUSION 



The semiconductor laser represents a good amplifier configuration which exhibits stable as well as strong chaotic 
turbulent states. The synchronization of semiconductor lasers in their chaotic regime is still a open area of research. 
The corresponding system of PDEs is, infinite dimensional and continuous coupling terms necessary for the syn- 
chronization. The coupling directly injects the output power from driving system to response system. Numerical 
simulations show that the complete synchronization is achieved although both the systems are chaotic in nature. 
We have used the dynamics of the synchronized states of the system, which play the roles of the transmitter and 
receiver sections, in colored image cryptography. The image can be securely transmitted and recovered over the 
secure communication channel. The process is applied to a colored image affecting the three color planes. The key 
space is vast enough that attackers may very hardly guess a key. Also, the speed of encryption is notable since the 
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key set preprocessed before hand. The second phase of encryption with a chaotic phase mask renders the system 
immune to common statistical attack, even if the initial conditions are identified thereby, forming a protective shield. 
Furthermore, the scheme can be extended by encrypting several image maps of the chaotic data set, thereby adding 
several layers of security. The proposed scheme has been tested exhaustively. Statistical tests re-instate that the 
scheme is also sensitive to its keys and it is almost impossible to infer any information about the original image if the 
chaotic set is tampered with. Overall, our method is robust and found to be resistant to most of the common threats. 
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